Draft version June 29, 2011 

Preprint typeset using I^T^X style cmulatcapj v. 11/10/09 



THE WHITE DWARF COMPANION OF A 2 M NEUTRON STAR 



o 

(N 

a 

(N 

P? 
43 

Oh- 
i < 

o ■ 



> 
OS 

m 

O 



13 



Varun B. Bhalerao and S. R. Kulkarni 

Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 

Draft version June 29, 2011 

ABSTRACT 

We report the optical discovery of the companion to the 2M Q millisecond pulsar PSR J1614-2230. 
The optical colors show that the 0.5 Mq companion is a 2.2 Gyr old He-CO white dwarf. We infer that 
M during the accretion phase is < 10 -2 M e dd- We show that the pulsar was born with a spin close to 
its current value, well below the rebirth line. The spin-down parameters, the mass of the pulsar, and 
the age of the system challenge the simple recycling model for the formation of millisecond pulsars. 
Subject headings: pulsars: individual: PSR J1614-2230 



1. PSRJ1614-2230 

PSR J1614-2230, a 3.15- ms pulsar, was discovered 
in a radio survey of unidentified EGRE T gamma ray 
sourc es using the Parkes Radio Telescope (|Hessels et al.l 
120051). Subsequently, X-ray emission from Newton XMM 
(jRoberts et al.l 120071 ) and 7-ray emission from Fermi 
Gam ma Ray Large Area Space Telescope ()Abdo et al.l 
was detected. Like most millisecond pulsars 
(MSPs), PSRJ1614-2230 is in a binary system. The 
circular orbit is consistent with the pulsar having un- 
dergone mass transfer and spun up. The mass function 
derived from pulsar timing indicated a companion with 
mass M 2 > 0.4 Mq (jHessels et al.ll2005H . 

The system recent ly came into prominence when 
iDemorest et al.l (|2010( ) reported the mass of the pulsar 
to be 1.97 ± 0.04 Mq. The detection of such a mas- 
sive neutron star (NS) places very strong constraints on 
the equation of state o f matter at extreme nuclea r dens i- 
ties (see, for example. lLattimer~ & Prakash 2004, 2005). 
The rather exquisite precision of this mass measurement 
was possible due to the orbit being almost perpendic- 
ular to the plane of the sky. As a result, the Shapiro 
delay caused by the companion is very large, result- 
ing in a precise estimate of the mass of the compan- 
ion, M 2 = 0.500 ± 0.006 Mq. The 8.7 day orbital pe- 
riod is significantly shorter than ~ 120 days expected for 
a low- mass X-ray binar y with such a massive secondary 
(|Rappaport et al.lll995D — suggesting a peculiar evolu- 
tionary history for this binary. 

Giv en the importance of the result of IDemorest et ail 
(2010) additional verification or consistency checks of 
physical parameters of PSR J1614-2230 can be expected 
to be of some value. A 0.5 Mq white dwarf (WD) at the 
inferred distance of PSR J1614-2230 (d ~ 1.2 kpc), even 
if a few Gyr old, is within the reach of present-day optical 
telescopes. It is this search for the WD that constitutes 
the principal focus of this Letter. 

2. OBSERVATIONS AT THE W. M. KECK OBSERVATORY 

We observed PSR J1614-2230 (Figured]) in g and R 
bands using the imaging mode of the Low Resolution 
Imagin g Spectrograph (LRIS) on the 10 m Keck-I tele- 
scope ( Oke et all Il995f ) , with upgraded recQ and blue 



Table 1 

Positions and Magnitudes of Reference Stars 



ID 




RA 


Declination 




m g 




a — 


16 h 14 m 


<5 - (-22°) 






A 


36 s 


98(8) 


29' 18". 78(11) 


17.94(8) 


18.73(3) 


B 


.36 s 


29(7) 


29'31".97(7) 


17.71(8) 


18.38(3) 


C 


35 s 


88(9) 


30' 18". 70(10) 


19.76(10) 


20.41(4) 


D 


34 s 


39(5) 


30' 12". 11(12) 


19.15(8) 


21.24(4) 


E 1 ' 


35 s 


66 


31'04".17 


21.38(9) 


22.04(5) 


F 


34 s 


92(16) 


30' 59". 46(2) 


20.18(9) 


21.36(5) 


G b 


37 s 


50 


30'43".69 


20.71(8) 


21.87(4) 


Q 


36 s 


47(7) 


30' 35". 90(4) 


(Saturated) 


17.20(3) 


R 


35 s 


92(2) 


30'30".45(1) 


20.10(8) 


20.94(4) 


S 


36 s 


50(10) 


30' 13". 93(13) 


(Saturated) 


17.75(3) 



http:/ /www2. keck. hawaii.edu/inst/lris/lris-red-upgrade- 
notcs.html 



Note. — Stars A-G were used as reference stars for photometry. 
Right ascension and declination were obtained from USNO-B1.0 
unless otherwise specified. 

a R band magnit udes calculated using SDSS magnitudes and 
IJester et al. (2005) transformation equations (Section 0. The 
numbers in parenthesis do not include a 0.03 mag uncertainty in 
absolute calibration due to the transformations. 
b Coordinates obtained from our final i?-band science image. The 
World Coordinate System for this image was calculated using a 
total of 33 USNO-B1.0 stars (Section |2j. 

cameras (IMcCarthv et alJl99llSteidel et al.l20 04L Sev- 
eral images were acquired at each target location, dither- 
ing the telescope by small amounts between each expo- 
sure. The observing conditions on UT 2010 May 15 were 
poor (seeing 1".4), so only data acquired on UT 2010 
July 8 (R band seeing 0".85 FWHM) were used in this 
analysis. The total exposure on this night was 960 s in 
the R band and 1010 s in the g band. The plate scale is 
0".135 pixel -1 for both cameras 

The images were processed using IRAF. After bias cor- 
rection and fl at fielding, cosmic rays were rejected using 
L . A . Cosmic rtvan Dokkurnl[200l . The images were then 
aligned with xregister and averaged to produce the fi- 
nal image for each band (see Figure [1]). A World Co- 
ordinate System was calculated using USNO-B1.0 stars 
in the field, with the imcoords package. Using 33 stars 
in the field, for the final i?-band image we obtained an 
rms error of 0".14 in right ascension (R.A.) and 0".20 in 
declination, adding up to a radial error of 0".24. For the 
final g-band image, the residuals were 0".13 and 0".18 for 
R.A. and declination, respectively, giving a total error of 
0".22. 
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Table 2 

Coordinates of PSR J1614-2230 at Different Epochs 



Method 


Epoch 


Ecliptic 
Longitude (A) 


Ecliptic 
Latitude (/?) 


RA 
a - 16 h 14 m 


Declination 
8 - (-22°30') 


Timing 
Timing 


J2005.63 
J2000.00 


245.78827556(5) 
245.78826025(12) 


-1.256744(2) 
-1.256697(5) 


36 s . 5051(1) 
36 s . 5034(2) 


31". 080(7) 
30". 904(19) 


Timing 

Chandra 3, 


J2007.32 
J2007.32 


245.78828015(6) 


-1.256758(2) 


36 s . 5056(1) 
36 s . 50(15) 


31". 132(9) 
31". 33(20) 


Timing 
LRIS 


J2010.51 
J2010.51 


245.78828886(11) 


-1.256785(5) 


36 s . 5067(2) 
36 s . 50(16) 


31". 23(2) 
31". 72(20) 



Note. — The proper motion of the pulsar as obtained from timing observations is fix = 9.79(7) masyr 1 , \xg = 
—30(3) mas yr —1 . This proper motion is used to estimate the ecliptic coordinates of the pulsar at the epoch of the Chan- 
dra and LRIS observations. The right ascension and declination are calculated from ecliptic coordinates using the Euler 
program in IDL. The equinox in all cases is J2000. 

a Bore-sight corrected coordinates. The target and the reference source "E" are detected in the Chandra image. We extract 
source coordinates using celldetect. We assume that the Chandra coordinate system and our i?-band coordinate solution are 
related by a simple offset with no rotation. Using the X-ray and i?,-band coordinates of source E, we calculate that the offset 
is ax — or = 0".21, &x — $H = — 0".20 and correct the target position using this offset. 



We measured fluxes with aperture photometry using 
the IDL APPH0T package. For each night, we measured 
the see ing (FWHM) and set the aperture to one seeing 
radius (Mighell 1999). We extracted the sky from an an- 
nulus 5-10 seeing radii wide. We had observed ne arby a 
Sloan Digital Sky Survey (SPSS: lYorFet"^II[2000l) fielcfl 
with the same settings as the science field. The cali- 
bration field was observed immediately after the science 
exposures, and had airmass 1.4 as compared to 1.6 for 
the target. We used the magnitudes of six stars from 
that field to calculate the photometric zero point and 
calibrated six reference stars in the science field (Figure 
[TJ Table [TJ. R band magnitudes were calculated using 
photo metric transformation^ prescribed bv lJester et al.1 
(2005) for stars with R c — I c < 1.15. The typical uncer- 
tainty in g-band magnitudes for stars with m g ~ 20 is 
0.03 mag. In R band we have Atur = 0.07 for tur ~ 20, 
including the uncertainty in the transformations. 

3. DETECTION OF AN OPTICAL COUNTERPART 

In the vicinity of the nominal pulsar position, we find 
a faint source (labeled "P" ) in the R band (Figure [I} . 
We do not detect anything within 1" of the target in the 
g band. The optical coordinates of this source, the tim- 
ing position of the pulsar and a proposed X-ray coun- 
terpart are summarized in Table [2] To compare this 
with the source location, we first have to correct for 
the 33 masyr -1 proper motion of the source. The LRIS 
source P is about 0".50 South of the pulsar position (ex- 
trapolated for the epoch of LRIS observations). Given 
the 0".24 (1-a) astrometric uncertainty of the optical im- 
ages, this position is consistent with the location of the 
pulsar. 

The density of objects brighter than star P in this field 
is 0.02arcsecond~ 2 . Using the 0".85 seeing FWHM as 
the mean diameter of stars, we calculate a false identifi- 
cation probability of 1%. The excellent astrometric co- 
incidence and the low probability of chance coincidence 
embolden us to suggest that star P is the optical coun- 
terpart of PSR J1614-2230. 

2 Calibration SDSS field: a = 17 h 19 m 10 s .10, 8 = -14°38'46".0. 

3 Photometric transforms: V = g — 0.59(g — r) — 0.01 ± 0.01 and 
V - R = 1.09(r - i) + 0.22 ± 0.03. 



Counterpart P is located only 4". 2 from the 16.3 mag 
star USNO-B1.0 0674-0429635, and is contaminated by 
the flux in the wings of its point-spread function. The 
proximity to this bright source will bias both the pho- 
tometry and the astrometry of the counterpart. To cal- 
culate the bias, we injected fake Gaussian sources with 
FWHM matched to the seeing and brightness compa- 
rable to the faint object. We then measured the coor- 
dinates and magnitude of the injected source using the 
same procedure as for the faint object. We find that for 
separations rj 4". 2, the recovered position is systemati- 
cally pulled towards the bright star by 0".l - 0".2. This 
is small enough that we do not apply this correction. The 
same exercise led us to derive the photometric bias. The 
de-biased i?-band magnitude of P is tur — 24.3 ±0.1. 

To infer properties of the WD, we need to calculate 
its absolute magnitude. To calculate the optical ex- 
tinction, we assume that the ratio Nn/N e is constant 
along the line of sight in the direction of PSRJ1614- 
2230. The dispe rsion measure (DM) f or PSR J1614-2230 
is 34.5 pc cm" 3 (|Demorest et al.ll2010Q, and the total D M 
in this direction is 104 pc cm" " 3 (jCordes fc Lazidl2002r i - 
about a factor of three higher than the pulsar DM. This 
implies that the pulsar is behind approximately one third 
of the galactic absorbing column. We can then scale the 
total R band extinct ion in this direction (Ar total = 0.65; 
ISchleeel et afl fl998l) to get A R = 0.22. We assume 
A ff = 0.47 /im for the g band3, and use the stan dard 
redde ning law with Ry = 3.1 to get A g = 0.34 ()Coxl 
1200th. 

IDemorest et al.l (|2010D use the DM to estimate that 
PSRJ1614-2230 is at a distance d = 1.2 kpc. For this 
distance, the extinction-corrected i?-band absolute mag- 
nitude is Mr = 13.7 ± 0.1. Furthermore, they place a 
lower limit of 900 pc on the distance, using pulsar paral- 
lax from timing measurements. The corresponding ab- 
solute magnitude is Mr = 14.3 ± 0.1. 

4. PULSAR AGE AND BIRTH SPIN PERIOD 

The measured mass (0.5 Mq) and the inferred ab- 
solute magnitu de (Mr 13.7) whe n applied to WD 
cooling models (jChabrier et al.l [20001) lead to an age of 

4 http:/ /www.sdss.org/dr6/instruments/imager/index.html#filters 
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Figure 1. Top: R band LRIS image of the PSR J1614-2230 field. The expected location of the target for epoch 2010.52 is shown by 
two horizontal lines just below the center of the image. The CCD shows considerable blooming in rows and columns around bright stars. 
The hexagonal mirror shape gives six diffraction spikes separated by 60° . To avoid contamination of the target by these artifacts from the 
bright star Q, we set the position angle to 300°. Bottom: i?-band (left) and <?-band (right) images of PSR J1614-2230. The target (P) is 
detected in the R band at a = 16 h 14 m 36 s .50, 8 = -22°30'31".72 and is marked with thick horizontal lines. There is a 0".33 offset between 
the expected and observed positions of the target (Table [2{. The target is not detected in the g-band, but the R band location is marked 
for reference. The diagonal streak in the R image is a bad CCD column. 
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Figure 2. Inferred parameters for PSR J1614— 2230 as a function of distance. The gray or dotted regions show the error bars on each 
parameter. Top left: The absolute ij-band magnitude of the WD. Top right: WD cooling age inferred from Chabrier ct al. (2000). The 
kink at 3Gyr may be a result of granularity of the tables. Bottom left: the intrinsic period derivative (P) in the pulsar frame, corrected 
for the Shklovskii effect. Since the pulsar cannot be spinning up, the values of P < are unphysical and are shown as the dotted area. 
The maximum distance to the pulsar is inferred to be 1540 pc. Lower right: the initial spin period of the pulsar. Birth periods slower than 
the current period are excluded. 
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twd ~ 2.2 Gyr. For such a WD, the expected absolute 
magnitude in other bands is Mb = 14.6 and My = 14.1. 
This gives m g ~ 25.0, consistent with our non detection. 

As per the standard evolutionary model for MSPs, the 
NS in PSR J1614-2230 was spun up by accretion from a 
low- or intermediate-mass companion star. The accretion 
stopped when the companion decoupled from its Roche 
lobe and became a WD. At this point, the MSP started 
spinning down by radiating energy. The spin-down age 
of the MSP is thus equal to the cooling down age of the 
WD. 

The period of the pulsar at birth (P ra it)> its spin-down 
age (r) , present-day period (P) and period derivative (P) 
are related to each other as follows: 



(n - 1)P 



P, 



P 



(1) 



where n is the "braking index" for the pulsar, with n = 3 
appropriate for a dipole radiating into vacuum. Thus, 
the period at birth is given by 



Pnit — P 



1 - 



tP(u - 1) 
P 



l/(r 



(2) 



The measured period derivative of a pulsar (P b s ) is al- 
ways higher than it s true period der ivative (P) owing to 
transverse motion ( ShklovskiH 1 1 9 70() . The corrected pe- 
riod derivative is P = P bs — Py, 2 d/c, where d is the dis- 
tance and \x is the proper motion. iDemorest et al.1 ()2010f ) 
measure \i = 32(3) masyr" 1 for PSR J1614-2230. Using 
the nominal distance d = 1.2 kpc, P = 4.8 x 10 _22 ss^ 1 . 

The DM inferred distance is quite uncertain, so it 
is useful to consider the dependence of all parameters 
on distance. Figure [2] shows the range of values for 
the P-band absolute magnitude (Mr), the WD cooling 
age (twd), the intrinsic period derivative in the pulsar's 
frame (P) and the initial spin period (Pmit). Since there 
is no energy injection to the pulsar now, it must be cur- 
rently spinning down. This implies P > and allows us 
to calculate an upper limit on the distance to the pul- 
sar: d max = 1540 pc (1-ct). The last panel in Figure [2] 
shows P init > 2.75 ms. We conclude that the pulsar in 
PSR J1614-2230 must have been born with a period close 
to the current observed value. 

The birth spin period of a NS is governed by an 
equilibrium between the ram pressure of the accret- 
ing material and the magnetic field. NSs spun up by 
accretion at the Eddington rate (M e dd) are reborn a s 
MSPs on the "rebirth line" (|Arzoumanian et al.lll999D . 
Figure [3] shows this rebirth line on a P-P diagram, 
along with curren t positions of various pu lsars from the 
ATNF database (|Manchester et all 12005). Also shown 
are rebirth lines for pulsars accreting at 10 _1 M c( jd and 
10~ 2 Af o dd- For a pulsar radiating as a dipole, the PP 
product remains constant through its lifetime. Thus, we 
can calculate the initial spin-down rate of the pulsar: 
Pnit = PP/Pnit, where P n it comes from Equation ([2]). 
If PSR J1614-2230 is at 900 pc, then it would have been 
born with P init = 3.1 x 10~ 21 ss _1 and P m t = 2.84 ms. 
This value is indicated with a star. The pulsar then 
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Figure 3. Pulsar P-P diagram. Period derivatives are corrected 
for the Shklovskii effect. Some interesting systems are labeled. 
Filled symbols are pulsars in globular clusters, whose P may be af- 
fected by the dense environment. Circles denote NS binaries with 
mass measurements. The dashed lines are the pulsar rebirth lines 
P = (M/M cdd )l.l X 10" 15 s- 4 / 3 P 4 / 3 HArzoumanian et al.lll999h 
for A//Af cdd = 1, KT 1 and 1CT 2 . The location of PSR J1614-2230 
on this plot depends on the distance (see Section[4jl. Were the pul- 
sar to be at 900 pc, it would be born at the star symbol and evolve 
toward its present-day location, the hollow circle marked "E". If 
instead it is at 1200 pc, the birth location is shown by the upward 
triangle, which also coincides with its present-day location. The 
solid curve passing through the star and the upward triangle de- 
notes all possible birth locations for PSR J1614— 2230. The system 
evolves through the shaded gray area to its present-day location 
on the vertical line through E and the upward triangle. 



evolves toward the lower right, to the circle marked "E" . 
Similarly, the birth parameters for the pulsar assuming 
d = 1.2 kpc are shown by the upward triangle — it is 
nearly coincident with the current parameters of the pul- 
sar for this distance. Other possible birth locations for 
the pulsar lie along the solid line passing through the 
star and triangle. The gray region shows all possible 
positions that PSR J1614-2230 can have occupied in its 
lifetime. It is clear that the pulsar was born well below 
the rebirth line, and the mean accretion rate during the 
spin-up phase (the final major accretion phase) was lower 

than 10~ 2 Af cd d- 

Two groups have run detailed simulations of t he evo- 
lution of PSR J1614-2230. The iLin et all ([2010^ model 
and the preferred model of lTauris et al.l (|2011[ ) are qual- 
itatively similar: the system begins as an intermediate 
mass X-ray binary consisting of a NS and a ~ 4M Q 
main-sequence secondary, which evolves to form a CO 
WD with an He envelope. The secondary accretes mass 
onto the NS in three phases. The first phase (Al) is 
a thermal timescale mass transfer at super-Eddington 
accretion rates, where the NS gains little mass. The 
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next phase (A2) is on a nuclear timescale (~ 35Myr), 
when the secondary is burning H in the core and en- 
velope. During this phase, the accretion rates are upto 
about a tenth of the Eddington limit. During the final 
accretion phase (phase AB), the secondary is burning 
He in its core and H in an envelope. This causes the 
radius of the donor star to expand, triggering accretion 
at near-Eddington rates for 5 - lOMyr. The NS gains 
the most mass during this phase. For typical NS param- 
eters, accreting 0.1 -0.2 is enough to spin them up 
millisecond periods (|Kiziltan &: Thor sctt 2010). Hence, 
the near-Eddington accretion in phase AB should spin 
the pulsar up all the way to the rebirth line. This is in- 
consistent with our inferred birth position of PSR J1614- 
2230. 

If the stellar evolution models are correct, then there 
is a pr oblem with the standard formation scenario of 
MSPs (|Radhakrishnan fc Srinivasanl 119821 : lAlpar et all 
1982). In essence, the birth period depends on fac- 
tors other than the magnetic field strength and accretion 
rates. For instance, iBildstenl ((19 98, 2003) proposes that 
accretion induces a quadrupole moment Q in the NS. 
The NS then loses angular momentum by gravitational 
wave radiation. Since the magnetic fields do not play 
any significant role in this model, the rebirth line be- 
comes irrelevant, and PSR J1614-2230 may start its life 
anywhere on the P-P diagram. In summary, the birth 
of PSR J1614-2230 away from the rebirth line requires 
reconsideration of angular momentum loss mechanisms. 

Going forward, on-going radio observations (pulsar 
timing and VLBI) should improve the parallax and 

thereby decrease the uncertainty in the inferred P and 
thus better locate the pulsar in the P-P diagram. The 
scattered light from light Q, while bothersome to the 
present observations, provide an opportunity to use 
adaptive optics to measure the near-IR fluxes of the WD. 
Grism spectroscopy with Hubble Space Telescope can po- 
tentially reveal the presum ed H+He layer posit ed by stel- 
lar evolutionary models of iTauris et al.l ([201 ll ). 
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